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Experimental study of the nuclear fusion reaction in charge-asymmetrical cfyi 3 He complex 
(d^ 3 He -> a (3.5 MeV) + p (14.64 MeV)) is presented. The 14.64 MeV protons were detected 
by three pairs of Si(dE — E) telescopes placed around the cryogenic target filled with the D2 + 3 He 
gas at 34 K. The 6.85 keV 7 rays emitted during the de-excitation of the d/i 3 He complex were 
detected by a germanium detector. The measurements were performed at two D2 + 3 He target den- 
sities, cp = 0.0585 and tp = 0.169 (relative to liquid hydrogen density) with an atomic concentration 
of 3 He 

C3u e — 0.0469. The values of the effective rate of nuclear fusion in d/i 3 He was obtained 
for the first time: X f = (4.5±§g) x 10 6 s -1 (<p = 0.0585); \ f = (6.9±|;g) x 10 5 s _1 (¥J = 0.168). 
The J = nuclear fusion rate in d/i 3 He was derived: A/ =0 = (9.7±|;£) x 10 5 s" 1 (<p = 0.0585); 
\j=° = (12.4±|;|) x 10 5 s" 1 (tp = 0.168). 

PACS numbers: 34.70.+C, 36.10.Dr, 39.10.+J, 82.30.Fi 
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I. INTRODUCTION 

The formation of muonic molecules of hydrogen iso- 
topes and their nuclear reactions have been the subject 
of many experimental and theoretical studies Q, 0, y, 0, 
IE IE Hi- As to the studies of formation of charge- 
asymmetrical muonic molecules like hfiZ (h = p,d,t, Z 
arc nuclei with a charge Z > 1) and their respective nu- 
clear fusion, the situation slightly different. What gave 
an impetus to study such systems was the theoretical 
prediction and experimental observation of the molecu- 
lar mechanism for charge exchange (MMCE) of p/i atoms 
on He nuclei |E . Essentially, the mechanism is re- 
duced to the following. Colliding with a He atom in a 
if-He mixture (H = H 2 , D 2 , T 2 and He = 3 Hc, 4 He), the 
muonic hydrogen atom forms a muonic complex hfiHc in 
the excited 2pa state. In the case of a deuterium-helium 
mixture, the complex may then decays from this state 
(see Fig. ^) via one of three channels 



dn + He [(dfiHc)*e-] + + e~ 
l 

[(^He)+e-]+7 (la) 
^ [(MHe)+e-]+d (lb) 

±* (/iHe)+ +d + e~. (lc) 
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If He = 3 Hc, fusion reactions may occur 

d/i 3 He a + n + p (14.64 MeV) (2a) 

^ /i 5 Li + 7 (16.4 MeV) . (2b) 

Thus, the fusion proceeds by the formation of a dfi atom, 
which, when incident on a 3 Hc atom, forms the c?^ 3 Hc 
molecular system. This molecule has two primary spin 
states, J = 1 and J = J&l; formation favors the former, 
fusion the latter I n Eqs. (jlal -c). A 7 is the (d^iHc)* 
molecular decay channel for the 6.85 keV 7-ray emission, 
A e for the Auger decay, and A p for the break-up process. 
The dfjMe molecule is formed with a rate A^He- The 
main fusion process, Eq. occurs with the rate A/, 

whereas the reaction (PJd), with the associated rate A/r 
has a branching ratio on the order of lO - ^ 4,5 - 1 [T^ . 

Interests in further study of charge-asymmetrical sys- 
tems was caused by first getting information on charac- 
teristics of the strong interaction in the region of ultralow 
energies. Secondly, it allows us to test the problem of 
three bodies interacting via the Coulomb law. More pre- 
cisely, these studies may allow us to 

- check fundamental symmetries and to measure the 
main characteristics of the strong interaction in 
the region of astrophysical particle collision ener- 
gies (~keV) in the entrance channel. It should be 
mentioned that nuclear fusion reactions in charge- 
asymmetrical muonic molecules are characterized 
by the same astrophysical range of energies ■ 

- test the calculation algorithm for rates of nuclear 
fusion reactions in /i-molecular complexes as well 
as for partial rates of decay of these asymmetrical 
complexes via various channels. 
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- solve some existing astrophysical problems. 

By now the experimental discovery of the MMCE has 
been confirmed in a number of experiments on study of 
muon transfer from h/j, to the He isotopes. 

Formation rates of the charge-asymmetrical duHe, and 

iljmiil and calculated plSSI^ l l^ lS i l il 
|32( with quite a good accuracy, and partial decay rates 
of such complexes were found. 



TABLE I: Experimental and calculated nuclear fusion rates, 
in s _1 , in the efyi 3 He complex. A/ is the effective rate of 



fusion reaction pkL \j 



and A 



are the rates of fusion 



reaction (j^Ji) in the <f/j 3 He complex in the J = and J = 1 



states, 


respectively. 




Refs. 


experiment 
[33] [34] 122 


[36] 


A/ 


< 7 x 10 7 < 1.6 x 10 5 < 6 x 10 4 


< 5 x 10 5 


Refs. 


theory 

[26] [37] [38] [39,40] [41J 


[42] 


\ .7=0 
X f 

\ J=1 
X f 


3 x 10 8 3.8 x 10 6 ~ 10 6 10 11 

to 6 


1.9 x 10 5 
6.5 x 10 2 



In the past five years interest in studying charge- 
asymmetrical complexes and in particular fusion in the 
cfyt 3 He system has revived. Tabic H] presents the calcu- 
lated fusion rates of deuterium and 3 He nuclei in the 
cfyx 3 He complex in its states with the orbital momenta 
J = and J = 1 and the experimental upper limits 
of the effective fusion rate, A/ , in the molecule averaged 
over the populations of fine-structure states of the <i/i 3 He 
complex. 

The experimental study of nuclear fusion in the d/i 3 Hc 
molecule is quite justified as far as detection of the pro- 
cess is concerned because there might exist an interme- 
diate resonant compound state 5 Li* leading to the ex- 
pected high fusion rate which results from quite a large 
value of the ^-factor for the ci 3 He reaction 43] . However, 
as follows from the calculations presented in Table Q] the 
theoretical predictions of the fusion rate in this molecule 
show a wide spread in value from ~ 10 5 s _1 to 10 11 s — 1 . 

The nuclear fusion rate in muonic molecules is usually 
calculated on the basis of Jackson's idea which allows 
the factorization of nuclear and molecular coordinates. In 
this case the nuclear fusion rate A„/ is given by 



A 



S 



nf 



(nMZ l Z 2 ) 



x |*. e (0)| 



(3) 



which is defined by the astrophysical (S-factor, the re- 
duced mass of the system M, the charges of nuclei in the 
muonic molecule Z\ and Z 2 , and the three-body system 
wave function ^ sc (0) averaged over the muon degrees of 
freedom and taken at distances comparable with the size 
of the nuclei, i.e., for r — > because of the short-range 
nature of the nuclear forces. 



It should be mentioned that, strictly speaking, asym- 
metrical muonic molecules [Z\ ^ Z 2 ) do not form bound 
states but correspond to resonant states of the continu- 
ous spectrum. In this case an analogue of Eq. © is given 
in Ref. Ull as 



A n ./ 



1 Mk ( 



(itMZ 1 Z 2 ) 21 + 1 4tt 



r|* sc (0)| 2 , (4) 



where I is the orbital quantum number of the resonant 
state, k is the relative momentum corresponding to the 
resonant energy, T is the width of the molecular state 
and ^ sc (0) is the wave function for the state of scatter- 
ing at resonant energy. In the limit of a very narrow 
resonance when r — > Eqs. © and (@J coincide. How- 
ever, one should take into account the asymptotic part 
of the wave function responsible for an in-flight fusion, 
including the possible interferences between the resonant 
and nonrcsonant channels. 

Let us briefly discuss the calculated nuclear fusion rates 
in the d/j, 3 He reaction presented in Table [I] The value 
given in Refs. H||[33| were given with some references to 
a calculation by Kamimura but without any references 
to the calculation method. In Ref. (3^ the author used 
a small variation basis and the experimental value of the 
astrophysical factor S « 6.32 MeV x b and found the 
nuclear fusion rate in the cfyx 3 He molecule in the ,7 = 
state to be 3.8 x 10 6 s" 1 . 

In Refs. [3^ . l4p| | the nuclear fusion rate in the d^ 3 Hc 
complex from the J = state was calculated by vari- 
ous methods. Since the nuclear fusion rate in the lsa 
states of the d/x 3 He molecule is much higher than the 
fusion rate form the 2pa state (because of a far smaller 
potential barrier), the undcr-barricr 2pa — > lsa transi- 
tion was calculated with finding the transition point in 
the complex r-planc. This procedure is not quite un- 
ambiguous and therefore the nuclear fusion rate in the 
cfyi 3 He molecule was calculated in an alternative way by 
reducing it to the S'-factor and using experimental data 
on low-energy scattering in 3 He (dp) 4 He reactions from 
Ref. |39l | . However, the procedure of an approximation 
of the experimental data for the ultralow energy region 
leads to some ambiguity of the results. The results of 
the calculation by the above two methods may differ by 
a factor of five for the i^ 3 He molecule and by a factor of 
three for the d/i 3 ile molecule in question |40j . 

The highest nuclear fusion rate was obtained in 
Ref. ^lj. Unlike the case in Ref. ^(j; where the bar- 
rier penetration factor in the 2pa — > lsa transition was 
evaluated, in Ref. the contribution from the lsa state 
to the total wave function for the at small intcrnuclcar 
distances r was determined. The determination of the 
contribution from this state to the total mesomolecule 
wave function at small distance requires the solution of 
a multichannel system of differential equations, which is 
a complicated problem because of the singularity of the 
expansion coefficients at small distances r — > 0. As to 
the results of Ref. ^3 given in the last column of Ta- 
ble U it is difficult to judge the calculation method used 
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because the method for calculation the wave function at 
small distances was not presented in the paper. 

Different results of calculations of the fusion rate in 
the <i/^ 3 Hc molecule reflect different approximations of 
the solution to the Schrodinger equation for three parti- 
cles with Coulomb interaction. The main uncertainty is 
associated with the results at small distances and hence 
follows the spread of the calculated values for the nuclear 
fusion rate in the d[i 3 He molecule given in Table Q] When 
the adiabatic expansion is used, the important problem 
of convergence of this expansion at small distances is usu- 
ally ignored. Such problems vanish if the direct solution 
of the Faddeev equations in the configuration space is 
performed in Refs. 0, EE 113- F° r this reason the cal- 
culation of the fusion rate in the d/x 3 He molecule using 
Faddeev equations in order to adjudge discrepancies be- 
tween different theoretical results becomes very actual 
problem. 

Much less has been done to study the nuclear fusion 
reaction in the cfyx 3 He experimentally. The estimations of 
the lower limit for the fusion reaction rate, has been 
done by a Gatchina - PSI collaboration using an ioniza- 
tion chamber [UIHIIl Their results ( see Tabled) dif- 
fer by several orders of magnitude. Another experiment 
aimed to measure the effective rate, A f tP , of nuclear fu- 
sion reaction (J2^,) was performed by our team [H. A 
preliminary result, also as estimation of lower limit, is 
shown in Table [I] 

The purpose of this work was to measure the effective 
rate, A/, of nuclear fusion reaction m the G?/i 3 He 
complex with the formation of a 14.64 MeV proton at 
two D2 + 3 Hc mixture density values. 



II. MEASUREMENT METHOD 

Figure n shows a slightly simplified version of the ki- 
netics to be considered, when negative muons stop in 
the D2 + 3 He mixture. The information on the fusion 
reaction rate in the efyx 3 He complex can be gained 
by measuring the time distribution, dN p /dt, and the to- 
tal yield, N p , of 14.64 MeV protons. These quantities 
are derived from the differential equations governing the 
evolution of the J = 1, states of the cfyt 3 He molecules 

Establishing the time dependence of the number of 
dn 3 Hc molecules, A^ 3He (f), for the two possible states 
J is sufficient to predict the time spectrum of the fusion 
products. In the following, we will include the effective 
transition rate Am of the d/i 3 Hc complex between the 
states J = 1 and J = 0. The Am transition is important 
if the Aj and A° rates differ strongly from one to another, 

and an appropriate value of Am permits the two rates to 
be measured. This possibility can be checked by measur- 
ing the fusion using different concentrations and densities 
which should also help clear up the questions surround- 
ing the mechanism of the Axo transition |48j . which is 
predicted to scale nonlinearly with the density. 

There is a direct transfer rate from ground state d/i's 



to 3 Hc's but that rate is about 200 times smaller than 
the A d 3 Ho rate and will be ignored [49|. No hyperfine 
dependence on the A^He formation rate is expected since 
the molecular formation involves an Auger electron and 
bound state energies of many tens of electron volts [ij. 
Using the expectation that the cfyx 3 He is formed almost 
exclusively in the J = 1 state, the solution for the fusion 
products from the J = and J = 1 states is relatively 
straightforward given the dfx population. The recycling 
of the muon after cfyx 3 He fusion will be ignored due to 
the extremely small probability of the fusion itself, and 
thus the system of equations decouples into the <i/i 3 He 
sector, and the dci-fusion sector (where cycling will be 
considered). Since there is no expectation of a J = 
to J = 1 transition, i.e., A01, the d[i 3 Hc sector is easily 
solved. 

Formation of dfid molecules from a dfi in hyperfine 
state F = 3/2 and F = 1/2 is given by the effective 
rate Xf, whereas the branching ratio f3p an d sticking 
probability Ud model the number of muons lost from the 
cycle by sticking. In both the initial condition on the 
number of d/j, atoms, and in the cycling efficiency after 
dd fusion, qi s represents the probability for a d[i atom 
formed in an excited state to reach the ground state • 
Finally, Wd, represents the probability that the muon 
will be captured by a deuterium atom given that there 
are both D? and 3 Hc in the mixture: 



W d = 



Xn 2 



c d + Acs He A r 



A'X. 



3 He 



(5) 



where Cd and C3 Hc are the deuterium and helium atomic 
concentration. A is the relative muon atomic capture 
probability by a 3 He atom compared to deuterium atom, 
and A' is the same ratio measured with respect to gas 
fraction concentrations (X). An previous experimental 
measure exists for D 2 + 3 He (A = 1.7 ± 0.2) [Sill EL 
and theoretical calculations for A' have been made by 
J. S. Cohen H^: for D 2 + 3 Hc : A' = 0.78 and 
for HD + 3 He : A' = 0.68. Our gas mixtures have 
C3 Hc = 0.0496(10) and thus A3 Hc = 0.0946(20). By 
atomic concentration, and using the experimental value, 
we get Wd = 0.92(2). Using theory and the gas fraction 
the result is the same, Wd = 0.92. Using our own exper- 
iment A = 1.67^0 33, to determine Wd leads also to 
the exact same value. 

The differential equations governing the evolution of 
the J = 1, spin states of the d/i 3 Hc molecules are (see 
Fig.HJ: 



dt 
dt 



= +tpc3 Hc \ d 3 Hc N dfl - A^A^hc (6) 



(7) 



where A^ is the number of dfj, atoms and with the defi- 
nition 



As = (Ao + A p 7=1 + A^ =1 + A e J=1 + Xj= x ) (8) 

> (9) 



(Ao + A p / = + Af + Af=° + A / 
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(1-W d ) 



((jid)u,Fj— 


qis 








(1- 


/ M 3 He J 


\ He 
■''cap 






[(d/x 3 He)*e _ ] + + e" 
2pa J=l 



I 



1— 



/j 3 He ^ +?i 
- 3 Hc + /J" + n (2.45 MeV) 



1-/8* 



t + +p(3.02 MeV) 



a + /j~ +p (14.64 MeV) 
/AA + 7(16.4 MeV) 



\J=l \J=l \J=1 



[/j 3 He + d] 



[(d/j 3 He)*e"]+ +e 
2po- J=0 



J 



\J=0 X j=o \J=0 



FIG. 1: Scheme of mu-atomic and mu-molecular processes occurring at stops of negative muons in the D2 + 3 He mixture 
and 



X d fi — A + ^C3 He A d 3 He 

+ ^c d A F [1 - W d q ls {l - (3 F uj d )} . (10) 

The yield for protons between two given times after the 
muon arrival, t\ and i 2 , is: 

Y p (h,t 2 ) = Y^(t 1 ,t 2 ) + Y°(t 1 ,t 2 ) 

AfD/Hc A/ <pc 3Kc \ d 3 Hc W d gi s e Y e p 
Ay; A dfl 



where the difference in time exponents has been defining 
as the yield efficiency: 



and with the effective fusion rate defined as 

A10 



A 



— A^ 




A10 + A° 
A10 + AM 
A10 + / 



+ A 



7=0 
/ 



A10 + A° 



(12) 

(13) 
(14) 



In the above equations, N^ Hc is the number of muons 
stopped in the D 2 + 3 He mixture and <p is the mixture 



atomic density relative to the liquid hydrogen density 
(LHD, N = 4.25 x 10 22 cm 3 ). 

When protons are detected in coincidences with muon 
decay electrons, later on called the del-e criterion, the 
fusion rate from Eq. I|ll|) takes the form 



X f = 



Y p (ti,t2)X dfl X^ 



N, 



D/Hc 



Wd qis <P C3 He A,j3He£p e e £t £y 



(15) 



where e e is the detection efficiency for muon decay elec- 
trons and St defined as 



-Ant, 



(16) 



is the time efficiency depending on the interval during 
which we accept the muon decay electrons. Note that 
Eqs. (|11H15|) arc valid when the proton detection times 
are t ^> 1/As • The values e p and As are found through 
calculation. Note an important feature of this experimen- 
tal setup: Xf is found by using the experimental values 
of Ad M , s e , W d , A d 3 Hc , and q ls . 

The information on these quantities corresponds to the 
conditions of a particular experiment and is extracted 
by the analysis of yields and time distributions of the 
6.85 keV 7 rays from reaction (|la|) . prompt and delayed 
x rays of /i 3 Hc atoms in the D 2 + 3 He mixture and muon 
decay electrons. The quantity A^He is determined from 
Eq. JTUJ) where F = 0.58, uj d = 0.122(3) are taken from 
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Eri 



10 cm 



FIG. 2: Apparatus used in the fiE4 area. The view is that 
of the incoming muon. Note that the Tl and TO scintillators 
are not shown. The labels are explained in the text. 



Refs. 0. X F = 0.05 x 10 6 s~ 1 is taken from Ref. 
The rate Xdp, is the slope of the time distribution of 7 ray 
from reaction i|la|) . 

The procedure of measuring qi s , A^hcj Wd , e e , A and 
A 7 (the partial probability for the radiative <i/i 3 He com- 
plex decay channel) as well as our results are described 
in detail in our previous work |56j |. 



III. EXPERIMENTAL SETUP 



decay electrons were detected by four pairs of scintilla- 
tors, E(/p, Edoi Eflj and E^s, placed around the vac- 
uum housing of the target. The total solid angle of the 
electron detectors was « 17%. The cryogenic target was 
located inside the vacuum housing. The design of the 
target is described in detail in Refs. |57ll58j. 

The analysis of the 6.85 keV 7-ray time distributions 
allows us to determine the disappearance rate, A^, for 
the d/j, atoms in the D2 + 3 He mixture. Note that the 
presence of a signal from the electron detectors during a 
certain time interval (the del-e criterion) whose beginning 
corresponds to the instant of time when the Ka 7 Kft, and 
K'f lines of [iHe atoms is detected makes it possible to 
determine uniquely the detection efficiency for muon de- 
cay electrons. When the del-e criterion is used in the 
analysis of events detected by the Si(dE — E) telescopes 
one obtains a suppression factor of 300 — 400 of the back- 
ground, which is quite enough to meet the requirements 
of the experiment on the study of nuclear fusion in the 
cfyr^Hc complex. 

Our experiment included two runs with the D 2 + 3 He 
mixture. The experimental conditions are listed in Ta- 
ble[nl In addition, we performed different measurements 
with pure D2, 3 He, and 4 He at different pressures and 
temperature. Details are given in Refs. [23j |. 



TABLE II: Experimental conditions for the D2+ He mixtures 
with an atomic concentration of helium C3 Ho — 0.0496. is 
the number of muon stopped in our apparatus. 



Run 




T 


P 


f 






[MeV/c] 


[K] 


[kPa] 


[LHD] 


[10 9 ] 


I 


34.5 


32.8 


513.0 


0.0585 


8.875 


II 


38.0 


34.5 


1224.4 


0.1680 


3.928 



The experimental layout (see Fig. [2J was described in 
details in Refs. j2jj|5|||53. The experimental facility was 
located at the //E4 beam line of the PSI meson factory 
(Switzerland) with the muon beam intensity around 2 x 
10 4 s~ 1 . After passing through a thin plastic entrance 
monitoring counter muons hit the target and stopped 
there initiating a sequence of processes shown in Fig. ^ 
The electronics are protected from muon pileup within a 
±10 /xs time gate so pileup causes a 30% reduction in the 
effective muon beam. Thus, we have a number of "good 
muons", called N^, stopping in our target. 

Three pairs of Si(dE — E) telescopes were installed 
directly behind 135 /im thick kapton windows and a 
0.17 cm 3 germanium detector behind a 55 /im thick kap- 
ton window to detect the 14.64 MeV protons from re- 
action (J2Jl) and the 6.85 keV 7 rays from reaction (|laf> . 
respectively. The Si telescopes with a 42 mm diameter 
were made of a 4 mm thick Si(-E) detector and a thin, 
360 /im thick, Si(dE) detector, respectively. An assembly 
of Si detectors like that give a good identification of pro- 
tons, deuterons, and electrons based on different energy 
losses of the above particles in those detectors. Muon 



The germanium detector was calibrated using 55 Fe and 
57 Co sources. The Si(dE — E) detectors were calibrated 
using a radioactive 222 Rn source. Before the cryogenic 
target was assembled, a surface saturation of the Si(<iE) 
and Si(E) detectors by radon was carried out. The 222 Rn 
decay with the emission of alpha-particles of energies 5.3, 
5.5, 6.0, and 7.7 MeV were directly detected by each 
of the Si detectors. The linearity of the spectrometric 
channels of the Si detectors in the region of detection 
of protons with energies 8 — 15 MeV was checked using 
exact-amplitude pulse generators. 



IV. ANALYSIS OF THE EXPERIMENTAL 
DATA 

A. Determination of the d/i 3 Hecomplex formation 

rate 

By way of example Fig. [3] shows energy spectra of 
events detected by the germanium detector in run I with- 
out and with the del-e criterion. The rather wide left 
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peak corresponds to the 7 rays with an average energy of 
6.85 keV and the three right peaks correspond to the Ka, 
Kfl, K-f lines of /iHe atoms with energies 8.17, 9.68, and 
10.2 keV, respectively. As seen in Fig. [3] the suppression 
factor for the background detected by the germanium 
detector with the del-e criterion is of the order of 10 3 . 

8000 



6000 



4000 



2000 



o 
u 



7 8 9 

Energy [keV] 



10 



1 1 



FIG. 3: Energy spectra of the events detected by the Ges 
detector in run I without (open circles) and with (full circles) 
the del-e criterion. 



Figure 0] shows time distributions of 6.85 keV 7 rays 
resulting from radiative de-excitation of the d/i 3 He com- 
plex in runs I and II. The distributions were measured 
in coincidences with delayed muon decay electrons. The 
experimental time distributions of 7 rays shown in Fig.0] 
were approximated by the following expression 



dN 7 



CPe 



-Ant 



L> 7 



(17) 



where B 7 , C 7 , and D 7 are the normalization constants. 
The second and third terms in Eq. (|17J) describe the con- 
tribution from the background. The analysis of the time 
distributions of the 6.85 keV 7 rays yielded values of Xdfi 
and thus the formation rates A^3 Hc . Results are given in 
Table ITTTI 

The systematic error is larger than the uncertainty of 
the result caused by various possible model of the back- 
ground, including the case where it is equal to zero (e.g., 
when time structure of the background is inaccurately 
known). We describe t he p rocedure of determining Arf3 Hc 
in more detail in Ref. 1391. 



TABLE III: Parameters used to determine the formation rates 
A d 3 He - The value Wd = 0.92(2) was used for both runs. 



Run 


51 s 










[10 6 s- 1 ] 


[10 6 s" 1 ] 


I 


0.882 (18) 


1.152(36)^(30)^ 


240(13) st at(15) SVBt 


II 


0.844 (20) 


2.496(58) siai (100) syst 


244(6) fltat (16) syst 



B. Number of muon stops in the D2 + 3 He mixture 

The number of muon stops in the D 2 + 3 He mixture 
was determined by analyzing time distributions of events 
detected by the four electron counters, We detailed this 
matter in Refs. [2^, I3TL l57j . Here it is pertinent to 
dwell upon some particular points in determination of 
this value. 

By way of example Fig. shows the time distribution 
of muon decay electrons measured in run I. To deter- 
mine the number of muon stops in the d/i 3 He the time 
distribution of the detected electrons, dN e /dt, is approx- 
imated by an expression which is superposition of four 
exponents and a background of accidental coincidences 

^ = A e M e- x " t +A e An e- XA - t +Al ic e- x ^ t +A e D e- Xat +B e , 

(18) 

where A e M , A\ u , A^ e and Afy, are the normalized ampli- 
tudes with 



AS 



KQiXoEe 



Al, Au, He, D, 



and 



Aai 
Aau 

A He 



Qai ■ Ao + Xclp 



— Qau • Ao + A, 
= A + A^ a e p , 



cap 



(19) 

(20) 
(21) 



are the muon disappearance rates in the different ele- 
ments (the rates are the inverse of the muon lifetimes 
in the target wall materials). In reality, Eq. I |18}l is 
an approximation of a more complex equation, which 
can be found in Ref. 153. The different rates are 



A = 0.455 x 10 6 s 



- 1 and XX 



2216(70) s" 1 The 



nuclear capture rates in aluminum and gold, A 
0.7054(13) x lOS" 1 and A^" = 13.07(28) x 10 6 s 



)6 C ,-1 

taken from Ref. I61l 
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FIG. 4: Time distributions of the 6.85 keV 7-quanta resulting 
from radiative de-excitation of the dfj?He complex obtained 
in coincidences with the muon decay electrons in run I. 
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take the new form 
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Time [us] 

FIG. 5: Time distributions of muon decay electrons measures 
in run I. The solid curves are the results of fitting its compo- 
nents (see Eq. 
background. 



23J): 1 - Au; 2 - Al; 3 - D 2 + J He; 4 - constant 



which take into account that muons are bound in the Is 
state of the respective nuclei when they decay. This fac- 
tor is negligible for helium but necessary for aluminum 
Q M = 0.993 and important for gold Q Au = 0.850 |6lj . 
The constant B e characterizes the random coincidence 
background. 

We denote iV M as the total number of muon stops in the 

D/He 



target, Nr- , iV„, and Np' as the numbers of muon 



rAl 

stops in Al, Au, and the gaseous D2 + 3 He mixture, re- 
spectively. Thus, we have the relation 



N n = N, 



Al 



N 



An 



/V D /Hc 



(22) 



Since the muon decay with emission electrons in the 
D2 + 3 Hc mixture take place from the Is state of the 
d/j, or /i 3 He atom, the third and fourth terms in Ea l|18(l 
will differ only by the values of the amplitudes A^ e and 
Afj because the slopes of both exponents are practically 
identical (A Ho = 0.457/^s -1 , A = 0.455 ^s" 1 )- In this 
connection the following simplified expression was used 
to approximate experimental time distributions of 



dt 



A 



- x A1 t 



Al' 



A 



Au c 



^D/Hc e ~ 



-Aau* _j_ 45 __ -A D/H c* _|_ _ge 

(23) 

Under our experimental conditions of runs I and II, we 
obtained the effective rates Xu/He = 0.4563 fis -1 and 
0.4567 /is" 1 , respectively. With these effective muon de- 
cay rates, the uncertainty in the calculated number of 
muon stops in the gaseous D2 + 3 Hc mixture is negligibly 
small as compared with the more rigorous calculation of 
this value by Eq. itT%|) . 

The amplitudes in Eq. (|19f) are expressed in terms of 
the factors clai, a Am and dD/Hei defined as the partial 
muon stopping in Al, Au, and D2 + 3 He mixture, 



N' 

tin 

N„ 



= 1 



Al, Au, D/He; (24) 



A\ = N^XoQtEedi 



(25) 



The electron detection efficiency, e e , of the detectors 
E[/p, Ebo, Eii/ and Ele was determined as a ratio be- 
tween the number of events belonging to the X-lines of 
the /x 3 He atoms, found from the analysis of the data with 
and without the del-e criterion, 



N x - 



(26) 



where N x _ e and N x are the numbers of events belong- 
ing to X-lines of the fi 3 He atom and detected by the 
germanium detector with and without coincidence with 
the electron detectors. The thus measured experimental 
value is electron detection efficiency averaged over the 
target volume. Table llVl presents the results. 

TABLE IV: Electron detection efficiencies, e e , in [%]. 



Run 






Detector 








E(7P 


Eh/ 


Edo 


Elb 


all 


I 


4.77(16) 


5.69(16) 


4.91(16) 


0.169(24) 


16.40(31) 


II 


4.53(15) 


5.89(18) 


4.88(14) 


0.114(39) 


16.34(39) 


£e 


4.65(12) 


5.79(12) 


4.89(12) 


0.148(23) 


16.37(22) 



The electron detection efficiency of the detector F,le 
is considerably lower than that of each of the other three 
electron detectors. This is because the material (Al, Fc) 
layer which the muon decay electron has to pass through 
in the direction of the detector E^£ is thicker than mate- 
rial layers in the direction of the other electron detectors. 

Table lists the values of the fraction of muons 
stopped in the D2 + 3 He mixture, ciD/Hei found from the 
analysis of the time distributions of the events detected 
by the four electron detectors in runs I and II. Note that 
when the au/He fraction, was calculated by Eqs. I|24|) 
and H25|) it was assumed that the electron detection effi- 
ciency by each of the detectors E[/p, Edo, ~Eri and E^b 
did not depend on the coordinates of the muon stop point 
in the target (be it in the target walls or in the D2 + 3 He 
mixture) . 



TABLE V: Fraction of muons stopped in the gaseous 
deuterium-helium mixture found by the absolute method. 
jyD/He j g numD er of muons stops in the D2 + 3 He gas 
mixture. 



Run 


a D/He 


N D,He 




[%] 


[10 9 ] 


I 


47.5(6) stot (30) s „ st 


4.216 


II 


66.6(10) st a*(39),„ 8 t 


2.616 



The systematic errors were determined as one half of 
the maximum spread between the a^/He values found 
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from analysis of the time distributions of the electrons 
detected by each of the electron detectors Ej/p, E^o, 
Ejii and File- Note that the fraction of rations stopped 
in gas, a D i He , is a result of simultaneously fitting all time 
distributions obtained with each of the electron detectors 
(and not a result of averaging all four distributions cor- 
responding to each of the four detectors). 

C. Determination of the detection efficiency for 
14.64 MeV protons 

To determine the proton detection efficiency, e p , of the 
three Si(dE — E) telescopes, one should know the distri- 
bution of muon stops over the target volume in runs I and 
II. The average muon beam momentum P M correspond- 
ing to the maximum fraction au /He of muons stopped 
in the D 2 + 3 He mixture in runs I and II was found by 
varying the muon beam momentum P M and analyzing the 
time distributions of the detected electrons by Eq. (|23[) . 
Next, knowing the average momentum P M and the beam 
momentum spread, we simulated the real distribution of 
muon stops in runs I and II by the Monte Carlo (MC) 
method [62J. The results of the simulation were used in 
another Monte Carlo program to calculate the detection 
efficiency of each pair of Si(dE — E) detectors for protons 
from reaction E3 ■ The algorithm of the calculation 
program included simulation of the muon stop points in 
the D2 + 3 Hc mixture and the d[i and ^ 3 He atom forma- 
tion points, the consideration of the entire chain of pro- 
cesses occurring in the D2 + 3 He mixture from the instant 
when the muon hits the target to the instant of possible 
production of 14.64 MeV protons in the fusion reaction in 
the cfyi 3 He complex. The calculation program took into 
account the proton energy loss in the gas target, kapton 
windows and Si(dE-E) detectors themselves (in the thin 
Si(dE) and thick Si(E) detectors). The proton detection 
efficiency s p was calculated at the gi s , W4, and values 
(see Table ITTTf) corresponding to our experimental condi- 
tions. The scattering cross sections of du atoms form D2 
molecules were taken from Refs. |64L I65L l66j . 

We ceased tracing the muon stopped in the target when 

a) the muon decays (n~ — * e~ + + v e ) 

b) the muon is transferred from the deuteron to the 
3 He nucleus with the formation of a 3 He/x atom 

c) nuclear fusion occurs in the c?/i 3 He complex 

d) a ddfi — > p + t + fi reaction proceeds in the ddfi 
molecule. 

Note that the algorithm of the program also involved the 
consideration of the background process resulting from 
successive occurrence of the reactions 

dfi + d^ddfi-> 3 Hc (0.8 MeV) + n 

+ 

d -> a + p(14.64MeV). (27) 



This reaction 1|27[) is called <i 3 Hc "fusion in flight" . 




0.0 0.2 0.4 0.6 0.8 1.0 

Energy of 3 He [MeV] 



FIG. 6: Cross section for the reaction 3 He + d — > 4 He + p in 
flight (reaction iJUJl) as a function of the 3 He nucleus-deuteron 
collision energy. The solid curve is the result of averaging the 
entire bulk of the presented experimental data. 

In our calculations we used the dependence of the cross 
section for reaction (|27() on the 3 He deuteron collision 
energy, averaged over the data of Refs. EE EE EE 
PFlI Fl ] . Figure displays the cross section dependence 
on the 3 He deuteron collision energy. The program also 
took into account the energy loss of 3 He nuclei in the 
D 2 + 3 He mixture caused by ionization of 3 He atoms and 
deuterium molecules. The time distributions of protons 
from reactions and (|27|1 under the same experimental 
conditions have completely different shapes in accordance 
with the kinetics of processes in the D2 + 3 He mixture. 

Figures. [7|and [8] show the calculated time dependen- 
cies of the expected yields of protons from reactions J5^) 
and (|27J) under the conditions of runs I and II. Thus, 
there arises a possibility of selecting a time interval of 
detection of events by the Si (dE — E) detectors where the 
ratio of the reaction J3^) and (|27() yields is the largest. 
This, in turn, makes it possible to suppress the detected 
background from reaction (|27|l to a level low enough to 
meet the requirement of the experiment on the study of 
nuclear fusion in the dfi 3 He complex. Tabic IVT1 presents 
the calculated values of some quantities describing kinet- 
ics of muonic processes in the D2 + 3 He mixture and the 
process of detecting protons from reactions (|2t) and $£7$ . 

TABLE VI: Calculated values of the quantities describing ki- 
netics of muonic processes in the D2 + 3 He mixture. The 
probabilities W3 Ua , W d/1 3n c , and W d 3 Ho are given per muon 
stop in the gaseous D2 + 3 He mixture. 



Run 


Ws He 




Wd3 He 




£ P 


JS 


Vp 


r\ s v s 




[io- 2 ] 


[io- 1 ] 


[IO" 5 ] 


[io- 1 ] 


[IO" 2 ] 


[IO" 2 ] 


[IO" 8 ] 


[IO" 8 


I 


2.60 


4.00 


2.735 


3.64 


3.40 


3.54 


2.26 


2.52 


II 


2.87 


5.16 


2.735 


2.06 


3.67 


3.47 


2.16 


2.72 



W3H0 is the total probability for the 3 He formation 
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FIG. 7: Calculated time distributions of the protons form reactions ^h) in runs I (left) and II (right). 
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FIG. 8: Calculated time distributions of the protons form reactions 1271 in runs I (left) and II (right). 



(i?3 He = 0.8 MeV) in the D2 + 3 He mixture, as a result 
of the fusion reaction in the dd\i molecule. W^He is 
the <i/x 3 He complex formation probability and W^Hc is 
the probability for d 3 He fusion in flight, following reac- 
tion 127|) . and W Me is the branching ratio of the muon 
decay via the pT — > e~ + + v e channel. e p and e'J 
are the detection efficiencies of one Si(dE — E) telescope 
for protons from reactions and i|^7|l . respectively. 
ijp and T}£* are the yields of protons from reactions J2t) 
and (|27|l detected by the Si(dE — E) telescope per muon 
stop in the gaseous D2 + 3 He mixture (the value of d/i 3 He 
fusion rate A/ = 10 6 s _1 was used for calculation of rj p ). 
There are some noteworthy intermediate results in the 
calculation of the detection efficiencies for protons from 
reactions and l|27|) . Table IVTT1 presents average en- 
ergy losses of protons on their passage through various 
material in the direction of the 8i(dE — E) detectors. 

Figures show the two-dimensional distributions of 
events detected by the Si(dE—E) detectors without coin- 
cidences with muon decay electrons in runs I and II. The 



x-axis represents the energy losses in the thin Si(dE) 
counters and the y-axis shows the total energies losses 
by the particle in both the Si(dE) and Si(-E) detectors 
connected in coincidence. The distributions of events in 
Figs. correspond to the detection of protons arising 
both from reactions (|2K) an d l|27|l and from the back- 



TABLE VII: Average energy losses, in MeV, of protons on 
their passage through various materials. 



Run 


D 2 + 3 He 


kapton 


Si(dE) 


Si(E) 




gas 


window 






I 


1.1 


0.6 


3.0 


10.1 


II 


3.5 


0.7 


3.7 


6.9 
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FIG. 9: Two-dimensional distributions of events detected by the Si(dE — E) telescopes in runs I (left) and II (right). 
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FIG. 10: Two-dimensional distributions of events detected by the Si(dE — E) telescopes in runs I (left) and II (right) with the 
del-e coincidence. 



ground reactions such as 

fx + 3 Hc — > p + 2n + Vfj, 
fi + Al -> Na* +p + n + v fl 
+P + v f _ l 

+p + 2n + (28) 
H + Fc -> Cr* +p + n + 

+p + 2n + v lx . 

In addition, the background which is not correlated with 
muon stops in the target (background of accidental coin- 
cidences) contributes to these distributions. 

Figures FTTH show the two-dimensional Si(dE—(E+dE)) 
distributions obtained in coincidences with muon decay 
electrons. As seen, the use of the del-e criterion leads 
to an appreciably reduction of the background, which 
in turn makes it possible to identify a rather weak ef- 



fect against the intensive background signal. To suppress 
muon decay electrons in the Si(dE — E) telescope, provi- 
sion was made in the electronic logic of the experiment to 
connect each of the electron detectors in anti-coincidence 
with the corresponding Si(dE — E) telescope. The choice 
of optimum criteria in the analysis of the data from the 
Si(dE — E) telescopes was reduced to the determination 
of the boundaries and widths of the time and energy in- 
tervals where the background is substantially suppressed 
in absolute value and the cffcct-to-background ratio is the 
best. To this end the two-dimensional Si(dE — (dE + E)) 
distributions corresponding to the detection of protons 
were simulated by MC method for runs I and II. On 
the basis of these distributions boundaries were deter- 
mined for the energy interval of protons from reaction 
(|2K) where the loss of the "useful" event statistics col- 
lected by the Si telescope would be insignificant. 

Figures ITU and IT21 show the two-dimensional Si(dE — 
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(dE + E)) distributions corresponding to the proton de- 
tection which were simulated by the MC method for runs 
I and II. Based on these distributions, we chose some 
particular proton energy intervals named AE^ when con- 
sidering the total energy deposited and 5E when looking 
only at the Si(dE) detector (see Table rvTTTf for further 
analysis. The regions of events corresponding to the in- 
tervals SE and AE^ are shown in the form of rectangles 
on the two- dimensional distributions presented on the 
Fig.HandEni 

It is noteworthy that the proton detection efficiencies 
given in Table IVTI correspond to these chosen proton en- 
ergy intervals for runs I and II. The next step in the data 
analysis was to choose a particular times interval of de- 
tection of events by the Si(dE — E) telescope. Figures HP 
and 1121 show the simulated time distributions of protons 
corresponding to the chosen energy loss intervals 5E, for 
the energy loss in the Si(cLE) detector and AEy, = E+SE 
the energy loss in both silicon detector. For the chosen 
proton energy intervals Table IVIIII presents the statis- 
tics suppression factors corresponding to different initial 
time,t t hr (with respect to the instant of the muon stop 
in the target) of the time intervals of detection of pro- 
ton events. These factors correspond to the ey value in 
Eq. JllJ). The data in Table lyTTTl are derived from time 
dependencies of the yields of protons from reactions 10}) 
and |(2ZJ) (see Figs. □ and EJ. 



According to the data given in Table IVTIIl we took the 
following time intervals At§\ (with isi the time for the Si 
signal to appear) for analyzing the events 

Ai Si (run I) : 0.7 < t Si < 2.2 (is 
At S i(runII): 0.4 < t Si <1.2^s. (29) 

Figures 1131 display the two-dimensional distributions of 
Si(dE — E) events obtained in coincidences with muon 
decay electrons in runs I and II with this time criteria 
imposed. With these time intervals Aisi and the proton 
energy loss AE% and SE intervals, the statistics collec- 
tion suppression factors for events from reactions 
and ijTT}) are 

kd^Hc = 2.9, k d 3 Hc = 11.2, Run I, 
fcd^Ho = 3.2, fc d 3 Hc = 12.1, Run II. 

(30) 

Another stage of the data analysis was the determina- 
tion of the number of events detected by the Si(dE — E) 
telescopes in runs I and II under the following criteria 

(i) the coincidence of signals from the Si telescopes 
and electron detectors in the time interval 0.2 < 
{t e — tsi) < 5.5 (is (t e is the time when the E de- 
tector signal appear). Such a requirement add the 
efficiency factor e t = 0.83 when determining the 
rates. 

(ii) the total energy release in the Si(dE) detector is 5E 
as given m Table IVTTT1 This particular SE interval 



will be called SE. For the thin and thick Si detector 
together, we choose the smallest interval, namely 
AEs = [11.7 - 14.2] MeV for run I and AE S = 
[8.0 - 13.4] MeV for run II. 

(iii) the time when the signal from the Si telescope ap- 
pears falls in the Aigj intervals. 



Table 11X1 presents the numbers of events N p detected in 
runs I and II under the above mentioned criteria. 

The contribution of the background events, N?f, given 
in Table Hxl from the reaction l|27|l is found in the follow- 
ing way. The expected number of detected protons from 
reaction (|27|) in runs I and II is calculated by 

f f _ N^a D/He W3 Hc W d m c £l f N si e e s t 

«d 3 Hc 

A^si is the number of 8i(dE — E) telescopes and l/fc d 3 He 
is the factor of background suppression by imposing the 
criteria (ii) and (iii). Using the values of a D / He and A M 
measured in runs I and II, the calculated values of Wa He , 
WdZHe, £pj Agi, fcjpHc, £ t, and Eq. i(3"TJ> . we obtained N£* , 
which is given in Table llXl Errors of the calculated 
arose from the inaccurate dependence of the cross sec- 
tions CT^He for the d 3 He reaction in flight on the 3 He 
deuteron collision energy and from the errors in the cal- 
culations of the detection efficiency of the Si telescopes 
for protons from reaction (|27fl . These errors were found 
by substitu ting various experimental <Jd?He{Ed?He) de- 
pendencies [ia S El Efl, El El- into the program 
for Monte Carlo calculation of the in-flight c? 3 He fusion 
probability W d 3 He . 

Now it is necessary to find the level of the accidental 
coincidence background by analyzing the experimental 
data from runs I and II. To this end the two-dimensional 
distribution of events detected by the Si(d£' — E) tele- 
scopes was divided into three regions which did not in- 
clude the separated region of events belonging to the 
process (j2t). Considering the boundaries of the inter- 
vals SE and AE^ of energy losses of the protons from 
reaction (|21}) we used three regions, A,B, and C, of the 
two-dimensional SE — AE-% distributions for determining 
the background level. The regions are given in Table Ixl 

The level N^ cc of the background of the accidental co- 
incidences of signals from the Si(dE — E) telescopes and 
the electron detectors for the given three region of the 
two-dimensional SE — AEs distributions and the corre- 
sponding suppression factor of the accidental background 
in the Si telescopes, rjsi-E, arc defined as 



i\jacc 



VSi-E 



N Si-E 



(32) 



(33) 



where Ag ; is the number of events detected by the three 
Si(dE — E) telescopes and belonging to the selected 
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FIG. 11: Two-dimensional distributions of Si(dE — (dE + E)) events obtained in run I by the Monte Carlo method and 
corresponding to detection of protons from reactions (left) and 12711 (right) in the time interval Atsi- 
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FIG. 12: Two-dimensional distributions of Si(dE — (dE + E)) events obtained in run II by the Monte Carlo method and 
corresponding to detection of protons from reactions (left) and 12711 (right) within the time interval Atsi- 



(SE — AEs) region of detection of protons from reaction 
(130- ^hi-E anci ^hi are t ne numbers of events detected 
by the ith Si(dE — E) telescope with and without del- 
e coincidences and belonging to the other (SE — AE^,) 
intervals. Note that the degree of suppression of the acci- 
dental coincidence background was determined not only 
by averaging the data obtained with the D 2 + 3 He mix- 
ture but also in additional experiments with the targets 
filled with pure 4 He, D 2 , and 3 He whose densities were 
if ~ 0.17, ip ~ 0.09, and <f ps 0.035, respectively. This 
guaranteed an identical ratio of stops in the target walls 
and in the gas in the experiments with 4 He, D 2 , and the 
D 2 + 3 He mixture {if — 0.168) and in the experiments 
with the D 2 + 3 Hc mixture (<p = 0.0585) and 3 He. Fig- 
ures^l^l and 1 161 display the two-dimensional distribu- 
tions of background events detected by the Si(dE — E) 



telescopes in the experiments with 4 He, D 2 , and 3 He. 

The values of rjsi-E an d N£ cc are given in Table llXl for 
runs I and II. The total numbers of detected background 
events, Np Ckg , which belongs to the analyzed region of 
energies (SE — AE-%) of protons from reaction (JSJi) and 
met the criteria (i)-(iii) were defined as 

N bckg = N ff + N acc ^ 

and are also given in Tablc llXl The uncertainties of Np Ckg 
include both statistical and systematical errors. 

Based on the measured values N v and the calculated 
values N^ ck s and following Refs. 0,[zl[zl, we found 
the yields of detected protons ,Y p , from reaction in 
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TABLE VIII: Time factors for reaction JSfc.) and reaction 12711 for the chosen intervals of energies of protons detected by the 
Si(dE — E) telescopes and detection beginning time, t±hr- 



Run 


A£ E 


5E 




Reaction J^Ji 


) 






Reaction 12711 






[MeV] 


[MeV] 






tthr, A* S 1 










tthr, /^S 1 








0.0 


0.2 


0.4 


0.7 


0.9 


0.0 


0.2 


0.4 


0.7 


0.9 


I 


[0 - oo] 


[O-oo] 


0.911 


0.684 


0.524 


0.350 


0.264 


0.989 


0.599 


0.388 


0.198 


0.131 




[11.7 - 14.2] 


[2.1 - 3.6] 


0.878 


0.659 


0.505 


0.337 


0.254 


0.438 


0.263 


0.171 


0.090 


0.058 


II 


[0 - oo] 


[O-oo] 


0.934 


0.543 


0.316 


0.129 


0.059 


0.996 


0.333 


0.114 


0.025 


0.009 




[8.0 - 13.4] 


[2.1-4.6] 


0.904 


0.525 


0.306 


0.125 


0.057 


0.752 


0.252 


0.084 


0.018 


0.006 
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FIG. 13: Two-dimensional distributions of events detected by the Si(dE — E) telescopes in runs I (left) and II (right) with the 
del-e coincidence and the time interval Atsi as define in Eq. 1291 . 



TABLE IX: Numbers of detected events, N v and iV/ / for the 
chosen 8E and AEt, intervals, taking into account the time 
intervals t e — tsi, and Atsi. Also the accidental coincidence 
background, Np CC , as well as the total background, Np Ckg . 



Run 


N p 


Ni f 




N acc 


jybckg 








[io- 4 ] 






I 


14 


3.8(2) 


4.2(9) 


2.5(5) 


6.3(6) 


II 


11 


2.4(1) 


2.4 (11) 


1.1(5) 


3.5(5) 



TABLE X: The three regions of division of two-dimensional 
(SE — AE-e) distributions. All energies are given in MeV. 





Region A 




Region 


B 


Region C 


Run 


AE S SE 




AE S 


5E 


AE S SE 


I 


0-11.7 3.6- 


6 


0-11.7 


-3.6 


14.2-25 1.8-6 


II 


0-8 4.6- 


6 


0-8 


-4.6 


13.6-25 1.5-6 



runs I and II. 



Yn 



7.7 
7.5 



+4.4 
-3.4 

-3.8 
-3.2 



run I 
run II 



(35) 

The errors of Y p are found in accordance with Refs. jx^, 
13 [HI dealing with analysis of small statistical samples. 
In view of Eq. l|ll|l and the measured values Y p , the ef- 
fective rate of nuclear fusion in the d/x 3 He complex is 
obtained from Eq. (|15|) . It can be written as 



V/- 



Y„ 



Nfj,a D / He Wdqi s <pc3 He \ d 3 Hc e p e e e t eY ' 



(36) 



The values of A/ and As corresponding to the conditions 
of runs I and II are given in Table IXII Using Eq. I|13|) 
and the measured effective rates of nuclear fusion and 
assuming that Aj- <C A° one can get hypothetical 

estimates of the partial fusion rate in the d[i 3 Ile complex 
in its states with the total orbital momentum J = 



A 



,/=o _ A/(Ai + A£) 



(37) 



Table IXll also presents the values for found in runs 
I and II. 
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FIG. 14: Two-dimensional distributions of events detected 
by the Si(dE — E) telescopes in a run with pure deuterium 
with the del-e coincidences and within the Atsi interval. The 
rectangle is the region corresponding to the energy intervals 
5E and AE^ for the run with the D2 + 3 He mixture at tp = 
0.168. 



FIG. 16: Two-dimensional distributions of events detected by 
the Si(dE — E) telescopes in the run with pure 3 He with the 
del-e coincidence and within the Aisi interval. The rectangle 
is region corresponding to the energy intervals 5E and AE^ 
for the run with the D2 + 3 He mixture at ip = 0.0585. 




TABLE XI: Effective rates of the 1 — > transition, A10, and 
the nuclear fusion rates in the d/^ 3 He complex. 



10 15 20 

Energy loss <5E— E 



60 



MeVl 



FIG. 15: Two-dimensional distributions of events detected 
by the Si(dE — E) telescopes in the pure 4 He with the del-e 
coincidences and within the Atsi interval. The rectangle is 
region corresponding to the energy intervals SE and AEs for 
the run with the D2 + 3 He mixture at tp — 0.168. 



The averages A° = 6 x 10 11 s" 1 and Xl = 7 x 10 11 s" 1 
(averaging over the data [H El El El HI El El) 

were used to get the values presented in Table IXII As 
to the effective rate for transition of the cfyx 3 He complex 
from the state with the angular momentum J = 1 to 
the state with J = 0, it was calculated with allowance 
for the entire complicated branched chain of processes 
accompanying and competing with the rotational 1 — 
transition (see Tabic |XT|) . The chain of these processes is 
considered in detail in Rcfs. 

EH El El El The effec- 
tive rates of nuclear fusion in the cfyt 3 He complex found 
by us in runs I and II coincide within the measurement 
errors. This is also true for the d 3 He fusion rates A^ =0 ob- 



Run 


A10 

[ io 1 ^- 1 ] 


\ J=0 

x f 

[ io^- 1 ] 


h 

[ lO'^s- 1 ] 


As 

[ io 1 ^- 1 ] 


I 

II 


5.2 
7.5 


q 7+5-7 
12.4±|-J 


4 -°-2.0 
6 Q+ 36 


6.54 
6.44 



tained by Eq. (|37|l . A comparison of the measured Aj =0 
with the theoretical calculations show rather goo d ag ree- 
ment with [3!|, a slight discrepancy with Refs. |38t It3| 
and considerable disagreement with Refs. |3ll l4l| . The 
cause of this disagreement is not clear yet as also is 
not clear the discrepancy between Aj =0 calculations in 
Refs. jH El El E3 f^ee Table D). Note that the the- 
oretical papers Refs. |37l l38l l39ll4lL yield estimates 
with a different degree of approximation. A correct com- 
parison of the experimental and theoretical A^ =0 is pos- 
sible only after carrying out some experiments with the 
D2 + 3 Hc mixture ruling out model dependence on the 
effective rate of transition of the d/x 3 He complex from the 
J = 1 state to the J = state. 

A comparison of the results of this paper with the ex- 
perimental results |35j reveals appreciable disagreement 
between them. The shortened form of presentation of the 
results j3f| does not allow us to find out sufficiently well 
the cause of this considerable disagreement. Note, how- 
ever, some results of the intermediate calculations which, 
to our mind, disagree with the real estimates of the cal- 
culated quantities. 

(1) According to Ref. J3{|, the fraction of the dfx atoms 
which were formed in the excited state under their 
experimental conditions and came to the ground 
state (per muon stop in the target) is = 0.8. 



15 



The quantity Cdp, is defined as 

C^ = \w P M^ + ^\ (38) 

where W p ,d is the probability for direct muon cap- 
ture by the HD molecule followed by formation of 
the muonic hydrogen atom or the excited dfi atom. 
qfs and are the probabilities for the transition 
of the p/i and d\i atoms from the excited state to the 
Is ground state. According to Refs. |2ll.l5ll80ll8ll|. 
under the Macv et. al experimental conditions the 
values of the quantities appearing in Eq. (|38[) were 
W p , d = 0.92 = 0.5, and = 0.8. Thus, as 
follows from our estimation, = 0.6 and not 0.8 
as stated. 

(2) The number of dp 3 He complexes formed in the 
course of data taking in their experiment was de- 
fined as 

N d ^ c = N,C d ,^f^, (39) 

and correspond to A^3 He = ( 4 -9 ± 0.4) x 10 s . 
According to our estimations, the quantities 
^dfj,3n c , Xdfi, Xpdfi (pdn molecule formation rate), 
and iVd M 3He had the values A^Hc = 1.32 x 10 6 s _1 
(<p = 0.0975 C3 Hc = 0.056, \% He = 2.42 x 

10 8 s" 1 ) m, 

ps 2.05 x 10 6 s" 1 (40) 

Xpdfj, = 5.6 x 10 6 s _1 , which yields A^„3 H c ~ 3.7 x 
10 s s" 1 instead of (4.9 ± 0.4) x 10 8 s^ f . 

(3) Their ionization chamber detection efficiency for 
protons from reaction was defined as e = esS T 
and found to be e = 0.082, where £5 = 0.13 is 
the selection factor for events detected in compli- 
ance with certain amplitude and geometrical cri- 
teria, e r = 0.63 is the time factor to take of the 
fact that the detected events were analyzed in the 
time interval 0.4 < t < 1.8 fxs. According to our 
estimation, e T = e _Ad ' jtl — e~ Xd ^ t2 = 0.44, because 
under their experimental conditions the dp, disap- 
pearance rate is Ad M ps 2.05 x 10 6 s -1 , t\ = 0.4 p,s, 
and tq. = 1.8 fis. 

As can be seen, taking into account only the above items 
alone the upper limit of A/ is, to our mind, apprecia- 
bly underestimated in the work of Maev et. al. An- 
other cause of this underestimation might be the im- 
proper background subtraction procedure because they 
determined the background level using information from 
earlier experiments |34j carried out under different condi- 
tions and at an experimental facility which was not com- 
pletely analogous. In addition, it is slightly surprising 



that the background from muon capture by 3 He nuclei 
with the formation of protons in the e nerg y re gion near 
14.64 MeV is estimated at zero in Ref. [35|(see p3|)- 

We believe that our A/ measurement results are reli- 
able, which is confirmed by stable observation of nuclear 
fusion in both runs with the D2 + 3 Hc mixture differing in 
density by a factor of about three. Nevertheless, as far 
as the experimental results obtained in this paper and 
in Ref. j35j are concerned, the things arc unfortunately 
uncertain and need clarifying. 

There is a point important for comparison of the cal- 
culated A^ =0 with the results of the previous experi- 
ments [3^ and this paper. Measurement of A^ =0 is indi- 
rect because it is determined by Eq. (|37|l with the calcu- 
lated effective rate for transition of the e?/i 3 He complex 
from the J = 1 state to the J = state. Therefore, A^ =0 

is not uniquely defined and greatly depends on A10, which 
in turn is determined by the chain of processes accom- 
panying and competing with the 1 — transition of the 
d/x 3 He complex. To rule out this lack of uniqueness in de- 
termination of Aj =0 and, in addition, to gain information 

on the effective 1—0 transition rate A10 and the nuclear 
fusion rate Xj =1 in the dfi 3 He complex in the J = 1 state, 
it is necessary, as proposed in Refs. [z3, 111, [Til , to carry 
out an experiment with the D 2 + 3 Hc mixture at least 
at three densities in the range (p = 0.03 — 0.2, where not 
only protons from reaction J2t) but also 6.85 keV 7 rays 
should be analyzed. Analysis of the results reported in 
this paper and in Ref. [35j makes it possible to put for- 
ward some already obvious proposals as to getting unam- 
biguous and precise information on important character- 
istics of /i-molecular (Xd^H e > A 10) and nuclear (A/, A^ =0 , 
A^ =1 ) processes occurring in the D 2 + 3 He mixture. It is 
necessary to conduct experiments at no less than three 
densities of the (HD + 3 He) or (H 2 + D 3 (l%) + 3 He) mix- 
ture with detection of both protons from reaction J2Ji) 
and 6.85 keV 7 rays, to increase at least three times the 
detection efficiency for protons e p and for muon decay 
electrons e e in comparison with the corresponding effi- 
ciencies in the present experiment. 
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